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Abstract 

Inclusive dijet and trijet production in deep inelastic ep scattering has been mea- 
sured for 10 < < 100 GeV^ and low Bjorken x, lO"'' < xbj < 10"^. The data 
were taken at the HERA ep collider with centre-of-mass energy ^/s = 318 GeV 
using the ZEUS detector and correspond to an integrated luminosity of 82 pb~^ 
Jets were identified in the hadronic centre-of-mass (HCM) frame using the fcr 
cluster algorithm in the longitudinally invariant inclusive mode. Measurements 
of dijet and trijet differential cross sections are presented as functions of Q^, xbj, 
jet transverse energy, and jet pseudorapidity. As a further examination of low- 
XBj dynamics, multi-differential cross sections as functions of the jet correlations 
in transverse momenta, azimuthal angles, and pseudorapidity are also presented. 
Calculations at 0{a^) generally describe the trijet data well and improve the 
description of the dijet data compared to the calculation at 
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1 Introduction 



Multijet production in deep inelastic ep scattering (DIS) at HERA has been used to 
test the predictions of perturbative QCD (pQCD) over a large range of negative four- 
momentum transfer squared, Q^, and to determine the strong coupling constant [1,2]. 
At leading order (LO) in as, dijet production in neutral current DIS proceeds via the 
boson-gluon-fusion {V*g qq with V = '-f, Z^) and QCD-Compton {V*q qg) processes. 
Events with three jets can be seen as dijet processes with an additional gluon radiation 
or with a gluon splitting into a quark-antiquark pair and are directly sensitive to 0{a1) 
QCD effects. The higher sensitivity to as and the large number of degrees of freedom 
of the trijet final state provide a good testing ground for the pQCD predictions. In 
particular, multijet production in DIS is an ideal environment for investigating different 
approaches to parton dynamics at low Bjorken-x, xbj [3]. An understanding of this regime 
is of particular relevance in view of the startup of the LHC, where many of the Standard 
Model processes such as the production of electroweak gauge bosons or the Higgs particle 
involve the collision of partons with a low fraction of the proton momentum. 

In the usual collinear QCD factorisation approach, the cross sections are obtained as 
the convolution of perturbative matrix elements and parton densities evolved according 
to the DGLAP evolution equations [4]. These equations resum to all orders the terms 
proportional to as In and the double logarithms In ■ In 1/ x, where x is the fraction 
of the proton momentum carried by a parton, which is equal to Xbj in the quark-parton 
model. In the DGLAP approach, the parton participating in the hard scattering is the 
result of a partonic cascade ordered in transverse momentum, p^- The partonic cascade 
starts from a low-p^ and high-x parton from the incoming proton and ends up, after 
consecutive branching, in the high-py and low-x parton entering in the hard scattering. 
This approximation has been tested extensively at HERA and was found to describe well 
the inclusive cross sections [5,6] and jet production [1,2,7,8]. At low xbj, where the phase 
space for parton emissions increases, terms proportional to In 1/x may become large 
and spoil the accuracy of the DGLAP approach. In this region the transverse momenta 
and angular correlations between partons produced in the hard scatter may be sensitive 
to effects beyond DGLAP dynamics. The information about cross sections, transverse 
energy, Et, and angular correlations between the two leading jets in multijet production 
therefore provides an important testing ground for studying the parton dynamics in the 
region of small Xbj. 

In this analysis, correlations for both azimuthal and polar angles, and correlations in 
jet transverse energy and momenta for dijet and trijet production in the hadronic {'~f*p) 
centre-of-mass (HCM) frame are measured with high statistical precision in the kinematic 
region restricted to 10 < < 100 GeV^ and 10~^ < a;Bj < 10~^. The results are compared 
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with pQCD calculations at next-to-leading order (NLO). A similar study of inclusive dijet 
production was performed by the HI collaboration [9]. 



2 Experimental set-up 

The data used in this analysis were collected during the 1998-2000 running period, when 
HERA operated with protons of energy Ep = 920 GeV and electrons or positron^ of 
energy E^. = 27.5 GeV, and correspond to an integrated luminosity of 81.7 ± 1.8 pb~^. A 
detailed description of the ZEUS detector can be found elsewhere [10, 11]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are measured in the central tracking detector (CTD) [12], which oper- 
ates in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consists of 72 cylindrical drift chamber layers, organised in nine superlayers covering the 
polar-angl^ region 15° < ^ < 164°. The transverse momentum resolution for full-length 
tracks can be parameterised as (j{pt)/pt = 0.0058pr © 0.0065 © 0.0014/pt, with pt in 
GeV. The tracking system was used to measure the interaction vertex with a typical 
resolution along (transverse to) the beam direction of 0.4 (0.1) cm and also to cross-check 
the energy scale of the calorimeter. 

The high-resolution uranium-scintillator calorimeter (CAL) [13] covers 99.7% of the total 
solid angle and consists of three parts: the forward (FCAL), the barrel (BCAL) and 
the rear (RCAL) calorimeters. Each part is subdivided transversely into towers and 
longitudinally into one electromagnetic section and either one (in RCAL) or two (in BCAL 
and FCAL) hadronic sections. The smallest subdivision of the calorimeter is called a 
cell. Under test-beam conditions, the CAL single-particle relative energy resolutions were 
a{E)/E = 0.18/ VE for electrons and (t{E)/E = 0.35/ VE for hadrons, with E in GeV. 

The luminosity was measured from the rate of the bremsstrahlung process ep — > e'-fp. The 
resulting small-angle energetic photons were measured by the luminosity monitor [14], a 
lead-scintillator calorimeter placed in the HERA tunnel a.t Z = —107 m. 



^ In the following, the term "electron" denotes generically both the electron (e ) and the positron (e+). 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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3 Kinematics and event selection 



A three-level trigger system was used to select events online [11, 15]. Neutral current DIS 
events were selected by requiring that a scattered electron candidate with an energy more 
than 4 GeV was measured in the GAL. The variable xbj, the inelasticity y, and were 
reconstructed offline using the electron (subscript e) [16] and Jacquet-Blondel (JB) [17] 
methods. For each event, the reconstruction of the hadronic final state was performed 
using a combination of track and GAL information, excluding the cells and the track 
associated with the scattered electron. The selected tracks and GAL clusters were treated 
as massless energy fiow objects (EFOs) [18]. 

The offline selection of DIS events was similar to that used in the previous ZEUS mea- 
surement [1] and was based on the following requirements: 

• i?g > 10 GeV, where E'^ is the scattered electron energy after correction for energy loss 
from the inactive material in the detector; 

• He < 0.6 and yjB > 0.1, to ensure a kinematic region with good reconstruction; 

• 40 < (5 < 60 GeV, where 5 = ~ Pz,i), where Ei and Pz,i are the energy and 
^-momentum of each final-state object. The lower cut removed background from 
photoproduction and events with large initial-state QED radiation, while the upper 
cut removed cosmic-ray background; 

• l^vtxl < 50 cm, where Z^tx is the Z position of the reconstructed primary vertex, to 
select events consistent with ep collisions. 

The kinematic range of the analysis is 

10 < g2 < 100 GeV^ 10"^ < XBj < 10-2 and 0.1 < y < 0.6. 

Jets were reconstructed using the kr cluster algorithm [19] in the longitudinally invariant 
inclusive mode [20]. The jet search was conducted in the HGM frame, which is equivalent 
to the Breit frame [21] apart from a longitudinal boost. 

The jet phase space is defined by selection cuts on the jet pseudorapidity, 77lab' 
laboratory frame and on the jet transverse energy, -E'j^hcm; HGM frame: 

-1.0 < r/S'^'^) < 2.5 and E^%^^ > 7 GeV, 4°|S > 5 GeV, 

where jetl,2(,3) refers to the two (three) jets with the highest transverse energy in the 
HGM frame for a given event. The dijet and trijet samples are inclusive in that they 
contain at least two or three jets passing the selection criteria, respectively. 
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4 Monte Carlo simulation 



Monte Carlo (MC) simulations were used to correct the data for detector effects, inef- 
ficiencies of the event selection and the jet reconstruction, as well as for QED effects. 
Neutral current DIS events were generated using the Ariadne 4.10 program [22] and 
the Lepto 6.5 program [23] interfaced to Heracles 4.5.2 [24] via Django 6.2.4 [25]. 
The Heracles program includes QED effects up to 0{a'^^). In the case of Ariadne, 
events were generated using the colour-dipole model [26], whereas for Lepto, the matrix- 
elements plus parton-shower model was used. The CTEQ5L parameterisations of the 
proton parton density functions (PDFs) [27] were used in the generation of DIS events 
for Ariadne, and the CTEQ4D PDFs [27] were used for Lepto. For hadronisation the 
Lund string model [28], as implemented in Jetset 7.4 [29,30] was used. 

The ZEUS detector response was simulated with a program based on Geant 3.13 [31]. 
The generated events were passed through the detector simulation, subjected to the same 
trigger requirements as the data, and processed by the same reconstruction and offline 
programs. 

The measured distributions of the global kinematic variables are well described by both 
the Ariadne and Lepto MC models after reweighting in [1] . The Lepto simulation 
gives a better overall description of the jet variables, but Ariadne provides a better 
description of dijets with small azimuthal separation. Therefore, for this analysis, the 
events generated with the Ariadne program were used to determine the acceptance 
corrections. The events generated with Lepto were used to estimate the uncertainty 
associated with the treatment of the parton shower. 

5 NLO QCD calculations 

The NLO calculations were carried out in the MS scheme for five massless quark fiavors 
with the program NLO JET [32]. The NLO JET program allows a computation of the 
dijet (trijet) production cross sections to next-to-leading order, i.e. including all terms up 
to 0{a1) iO{al)). In certain regions of the jet phase space, where the two hardest jets 
are not balanced in transverse momentum, NLOjet can be used to calculate the cross 
sections for dijet production at 0{al). It was checked that the LO and NLO calculations 
from NLOjet agree with those of Disent [33] at the 1-2% level for the dijet cross 
sections [34,35]. 

For comparison with the data, the CTEQ6M [36] PDFs were used, and the renormalisation 
and factorisation scales were both chosen to be {E^yicm~'<~Q'^)/^: where for dijets (trijets) 
-E'r.HCM is the average -Et.hcm of the two (three) highest -Et,hcm jets in a given event. 
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The choice of renormahsation scale matches that used in the previous ZEUS multijet 
analysis [1] . The strong coupling constant was set to the value used for the CTEQ6 PDFs, 
as{Mz) = 0.118, and evolved according to the two-loop solution of the renormahsation 
group equation. 

The NLO QCD predictions were corrected for hadronisation effects using a bin-by-bin 
procedure. Hadronisation correction factors were defined for each bin as the ratio of the 
hadron- to parton-level cross sections and were calculated using the Lepto MC program, 
which, at the parton level, gives a better agreement with NLOjet than Ariadne. The 
correction factors Chad were typically in the range 0.8 — 0.9 for most of the phase space. 

The theoretical uncertainty was estimated by varying the renormahsation scale up and 
down by a factor of two. The uncertainties in the proton PDFs were estimated in the pre- 
vious ZEUS multijets analysis [1] by repeating NLOjet calculations using 40 additional 
sets from CTEQ6M, which resulted in a 2.5% contribution to the theoretical uncertainty 
and was therefore neglected. 

6 Acceptance Corrections 

The Ariadne MC was used to correct the data for detector effects. The jet transverse 
energies were corrected for energy losses from inactive material in the detector. Typical 
jet energy correction factors were 1 — 1.2, depending on the transverse energy of the 
detector-level jet and the jet pseudorapidity. 

The measured cross sections were corrected to the hadron level using a bin-by-bin proce- 
dure. These corrections account for trigger efficiency, acceptance, and migration. Typical 
efficiencies and purities were about 50% for the differential cross sections, with correc- 
tion factors typically between 1 and 1.5. For the double-differential cross sections, the 
efficiencies and purities were typically 20 — 50%, with correction factors between 1 and 2. 

The cross sections were corrected to the QED Born level by applying an additional cor- 
rection obtained from a special sample of the Lepto MC with the radiative QED effects 
turned off. The QED radiative effects were typically 2 — 4%. 

7 Systematic uncertainties 

A detailed study of the sources contributing to the systematic uncertainties of the measure- 
ments has been performed. The main sources contributing to the systematic uncertainties 
are listed below: 
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• the data were corrected using Lepto instead of Ariadne; 

• the jet energies in the data were scaled up and down by 3% for jets with transverse en- 
ergy less than 10 GeV and 1% for jets with transverse energy above 10 GeV, according 
to the estimated jet energy scale uncertainty [37]; 

• the cut on -E^^hcm each jet was raised and lowered by 1 GeV, corresponding to the 
Et resolution; 

• the upper and lower cuts on ?7lab^*'''^'' were each changed by ±0.1, corresponding to the 
7] resolution; 

• the uncertainties due to the selection cuts was estimated by varying the cuts within 
the resolution of each variable. 

The largest systematic uncertainties came from the uncertainty of the jet energy scale, 
which produced a systematic uncertainty of 5 — 10%. For the trijet sample, altering the 
cut on -Et^hcm ^Iso produced a systematic uncertainty of 5 — 10%. The other significant 
systematic uncertainty arose from the choice of Lepto instead of Ariadne for correcting 
detector effects. This systematic uncertainty was also typically 5 — 10%. The other 
systematic uncertainties were smaller than or similar to the statistical uncertainties. 

The systematic uncertainties not associated with the absolute energy scale of the jets were 
added in quadrature to the statistical uncertainties and are shown as error bars in the 
figures. The uncertainty due to the absolute energy scale of the jets is shown separately 
as a shaded band in each figure, due to the large bin-to-bin correlation. In addition, there 
is an overall normalisation uncertainty of 2.2% from the luminosity determination, which 
is not included in the figures. 



8 Results 

8.1 Single-differential cross sections da/dQ'^, da/dx^j and trijet 
to dijet cross section ratios 

The single-differential cross-sections da/dQ"^ and da/dx^j for dijet and trijet production 
are presented in Figs. [U^a) and (c), and Tables [1]- HI The ratio (Ttrijct/o'dijct of the trijet 
cross section to the dijet cross section, as a function of and of xbj are presented 
in Figs. [U^b) and [I](d), respectively. The ratio o'trijet/o'dijct is almost independent, 
as shown in Fig. [TJ^b), and falls steeply with increasing xbj, as shown in Fig. [U^d). In 
the cross-section ratios, the experimental and theoretical uncertainties partially cancel, 
providing a possibility to test the pQCD calculations more precisely than can be done 
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with the individual cross sections. Both the cross sections and the cross-section ratios are 
well described by the NLOjet calculations. 

8.2 Transverse energy and pseudorapidity dependencies of cross 
sections 

The single-differential cross-sections da / dE^r^\^(jy^ for two (three) jet events are presented 
in Fig. [2j The measured cross sections are well described by the NLOjet calculations 
over the whole range in E^r^^^^^^ considered. 

The single-differential cross sections da/driQ^^ for dijet and trijet production are presented 
in Figs. [3]|^a) and[3]|^c). For this figure, the two (three) jets with highest -E^hcm were 
ordered in //lab- Also shown are the measurements of the single-differential cross-sections 
da / d\A'rf-^^^\, where |A?7j^^^j^| is the absolute difference in pseudorapidity of the two jets 
with highest E^r^'^^iCM (^^^ Figs. |3](b) and[3l^d)). The NLOjet predictions describe the 
measurements well. 

8.3 Jet transverse energy and momentum correlations 

Correlations in transverse energy of the jets have been investigated by measuring the 
double-differential cross-sections d^a / dx-Q^d^E^^^^^y^, where ^E^^^^^y^ is the difference 
in transverse energy between the two jets with the highest E^,^'^y{cm- "^^^ measurement 
was performed in xbj bins, which are defined in Table for dijet and trijet production. 
Figures m and [5] show the cross-sections d'^a / dx-B]dAE^^^^^y^ for all bins in x-q] for the dijet 
and trijet samples, respectively. 

The NLOjet calculations at 0{al) do not describe the high-Ai?;j!:*j^^j^ tail of the dijet 
sample at low xbj, where the calculations fall below the data. Since these calculations give 
the lowest-order non-trivial contribution to the cross section in the region ^E}^^^^^^ > 0, 
they are affected by large uncertainties from the higher-order terms in a^. A higher-order 
calculation for the dijet sample is possible with NLOjet if the region AE^^^^^y^ near zero 
is avoided. NLOjet calculations at 0{a^) for the dijet sample have been obtained for 
the region AE^^^^^y^ > 4 GeV and are compared to the data in Fig. HI With the inclusion 
of the next term in the perturbative series in a^, the NLOjet calculations describe the 
data within the theoretical uncertainties. The NLOjet calculations at 0{al) for trijet 
production are consistent with the measurements. 

As a refinement to the studies of the correlations between the transverse energies of the 
jets, further correlations of the jet tranverse momenta have been investigated. The correla- 
tions in jet transverse momenta were examined by measuring two sets of double-differential 
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cross sections: d'^a/dxB]d\Tip^^^(^\ and d'^a/dxB]d{\Apr^^^(^\/ (Q.E^^^^q-^)). The variable 
I^PthcmI is transverse component of the vector sum of the jet momenta of the two 
jets with the highest -Ej^^hcm- -^^^ events with only two jets jSp^ilj^Q^I = 0, and ad- 
ditional QCD radiation increases this value. The variable | Ap^"j^qmI/(2-E'thcm) 
magnitude of the vector difference of the transverse momenta of the two jets with the 
highest -E^HCM scaled by twice the transverse energy of the hardest jet. For events with 
only two jets \Ap^^^^\/ {2E^r^^^^^) = 1, and additional QCD radiation decreases this 
value. Figures E] - [9] show the cross-sections d'^a/dxBid\'E'pr^'^^(^\ and the cross-sections 
d'^a /dxBjd\Apr^^l^(^\/{2E^r^^^^j^) in bins of xbj for the dijet and trijet samples. 

At low XBj, the NLOjet calculations at O(a^) underestimate the dijet cross sections at 
high values of and low values of | Ap^"j^^^|/(2£';|?*j^(^p^). The description of the 

data by the NLOjet calculations at improves at higher values of x-Qy A higher- 

order calculation with NLOjet at 0{al) for the dijet sample has been obtained for the 
region |Ep^*j^j3^| > 4GeV, which is compared to the data in Fig. [6l and for the region 
I^PthcmI/(2-E'thcm) < 0.85, which is compared to the data in Fig. [HI With the inclusion 
of the next term in the perturbative series in a.^, the NLOjet calculations describe the 
data well. The NLOjet calculations at 0{al) for trijet production are consistent with 
the measurements. 



8.4 Azimuthal distributions of the jets 

Measurements of the double-differential cross-section d'^a /dxBid\A(j)^^(^^\, where |A</)i^^^^| 
is the azimuthal separation of the two jets with the largest -E-rHCM' dijet and trijet 
production are shown in Figs. [10] and [TT] for all bins in xbj. For both dijet and trijet 
production the cross section falls with |A0i^^^^|. The NLOjet calculations at for 
dijet production decrease more rapidly with |A0i^^^^| than the data and the calculations 
disagree with the data at low |A0i^^^^|. A higher-order NLOjet calculation at 0{ag) 
for the dijet sample has been obtained for the region |A0i^|^]^| < Sir/A and describes the 
data well. The measurements for trijet production are reasonably well described by the 
NLOjet calculations at 

A further investigation has been performed by measuring the cross-section d'^a / dQ'^dxB] 
for dijet (trijet) events with |A0i^^^^| < 27r/3 as a function of xb]- For the two-jet final 
states, the presence of two leading jets with |A0i^^^^| < 27r/3 can indicate another high-i?r 
jet or set of high-i^j- jets outside the measured rj range. These cross sections are presented 
in Fig. [121 The NLOjet calculations at 0{a\) for dijet production underestimate the 
data, the difference increasing towards low xbj. The NLOjet calculations at 0(a;|) are 
up to about one order of magnitude larger than the C(a|) calculations and are consistent 
with the data, demonstrating the importance of the higher-order terms in the description 
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of the data especially at low xbj- The NLOjet calculations at 0{a^) describe the trijet 
data within the renormalisation-scale uncertainties. 

9 Summary 

Dijet and trijet production in deep inelastic ep scattering has been measured in the phase 
space region 10 < < 100 GeV^ and 10~^ < xbj < 10~^ using an integrated luminosity 
of 82 pb~^ collected by the ZEUS experiment. The high statistics have made possible 
detailed studies of multijet production at low xbj. The dependence of dijet and trijet pro- 
duction on the kinematic variables and xbj and on the jet variables -E^hcm ^^"^ '7lab 
is well described by perturbative QCD calculations which include NLO corrections. To 
investigate possible deviations with respect to the collinear factorisation approximation 
used in the standard pQCD approach, measurements of the correlations between the two 
jets with highest E^T\iCM have been made. At low xbj, measurements of dijet produc- 
tion with low azimuthal separation are reproduced by the perturbative QCD calculations 
provided that higher-order terms ((9(a^)) are accounted for. Such terms increase the 
predictions of pQCD calculations by up to one order of magnitude when the two jets with 
the highest -E^hcm balanced in transverse momentum. This demonstrates the 

importance of higher-order corrections in the Iow-Xbj region. 
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( GeV') 


7^72 ^stat ^syst ^ES 

(pb/GeV^) (pb/GeV^) (pb/GeV^) (pb/GeV^) 


Cqed 


Chad 


10 - 15 
15 - 20 
20 - 30 
30 - 50 
50 - 100 


66.0 0.8 lil til 
41.4 0.6 til III 
26.2 0.3 til til 

14.0 0.1 1°:^ t\l 

5.82 0.06 


0.984 
0.968 
0.965 
0.955 
0.952 


0.866 
0.870 
0.876 
0.884 
0.887 



Table 1: The inclusive dijet cross sections as functions of . Included are the 
statistical, systematic, and jet energy scale uncertainties in columns 3, 4, ond 5, re- 
spectively. Column 6 shows the correction factor from QED radiative effects applied 
to the measured cross sections, and column 7 shows the hadronization correction 
applied to the NLOjet calculations shown in the figures. 



XBj X 10 


^stat ^syst ^ES 

(pb, xlO"^) (pb, xlO-^) (pb, xlO^^) (pb, xlO^^) 


Cqed 


Chad 


1.7-3.0 
3.0 - 5.0 
5.0 - 10.0 
10.0 - 25.0 
25.0 - 100.0 


85.3 1.7 til til 
113.8 1.5 til til 
83.1 0.8 til tri 
29 5 3 +°-^ +2-^ 

Z,J.O U.O _Q g _2 2 
Z.Oi U.UO _o.o7 -0.17 


0.987 
0.975 
0.969 
0.958 
0.948 


0.910 
0.887 
0.876 
0.876 
0.862 



Table 2: The inclusive dijet cross sections as functions of xb]. Other details as 
in the caption to Table 1. 
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( GeV') 


-J^ ^stat ^syst ^ES 

(pb/GeV^) (pb/GeV^) (pb/GeV^) (pb/ GeV^) 


Cqed 


Chad 


10 - 15 
15 - 20 
20 - 30 
30 - 50 
50 - 100 


'•^ '^■^ -1.3 -1.0 
^•^•J U.l/ _Qgg _0 52 

q 1 Q nil +0-27 +0.38 
'^■^^ ^-^^ -0.37 -0.38 

i.UO U.UU _o.ii „o.l9 

719 024 +°-°^^ +°-°^'^ 

U.(1J U.U^'l -0.027 -0.070 


0.991 
0.946 
0.969 
0.949 
0.956 


0.759 
0.776 
0.786 
0.794 
0.795 



Table 3: The inclusive trijet cross sections as functions of . Other details as 
in the caption to Table 1. 



xbj X 10 


^stat ^syst SeS 

(pb, xlO"^) (pb, xlO"^) (pb, xlO"^) (pb, xlO"^) 


Cqed 


Chad 


1.7-3.0 
3.0 - 5.0 
5.0 - 10.0 
10.0 - 25.0 
25.0 - 100.0 


14.7 0.7 t\i 
15.9 0.5 ^2.0 +1.9 

9.6 0.3 1°| 

OOP- Q1Q +0.21 +0.40 
O.Od U.iU _QiQ _Q37 

niQ9 nnT? +0.023 

U.iJZ U.UiO -0.020 -0.022 


1.00 
0.968 
0.961 
0.954 
0.95 


0.811 

0.796 
0.780 
0.785 
0.739 



Table 4: The inclusive trijet cross sections as functions of Xbj ■ Other details as 
in the caption to Table 1. 
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Variable 


Bin 


Boundaries 


A r;jctl,2 


2 


w t: vjrc V 

4-10 GeV 




3 


10 - 18 GeV 




4 


1 8 _ 1 no OeV 


1 ^Pt,hcm I 


T 

X 

2 


n — 4 Qpv 

4-10 GeV 




3 


10 - 16 GeV 




4 


16-100 GeV 




1 
2 


0-0.5 
0.5-0.7 




3 


0.7-0.85 




4 


0.85-1 


I^0hcmI 


1 
2 
3 
4 


- 7r/4 
7r/4-7r/2 
7r/2 - 37r/4 
37r/4 - TT 



Table 5: The bin edges used for the measurements of the jet correlations presented. 
For the trijet sample, the first two bins in |A0i^^^| are combined. 
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• ■ ZEUS 82 pb ' 

- NLOjet:0(ap®C,,,; 
NL0jet:0(a3)®C,^^: 




■ jet energy scalg uncertainty 

□ 1/16<|j2/(Q'+Et)<1 (^a) 



-t 1 1 1 1 — I — I — h 



H 1 1 1 1 1 1 — H 



c 0.15 



D 



50.14 



0.13 
0.12 
0.11 
0.1 

0.09 

0.2 



S 



-0.2 
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Q2 (GeV^) 



• ZEUS 82 pb ' 

tri-/dijets 
- NLOjet: 0(a^)/0(a2) 




jet energy scajg uncertainty: 
□ 1/16<h2/(q2+e^)<1 (^b)^ 



H 1 1 1 1 1 1 — I- 



20 30 40 50 10' 

Q2 (GeV^) 




Figure 1: Inclusive dijet and trijet cross sections as functions of (a) and (c) 
XBj- Figures (h) and (d) show the ratios of the trijet to dijet cross sections. The 
bin-averaged differential cross sections are plotted at the bin centers. The inner 
error bars represent the statistical uncertainties. The outer error bars represent the 
quadratic sum of statistical and systematic uncertainties not associated with the jet 
energy scale. The shaded band indicates the jet energy scale uncertainty. The pre- 
dictions of perturbative QCD at NLO, corrected for hadronisation effects and using 
the CTEQ6 parameterisations of the proton PDFs, are compared to data. The lower 
parts of the plots show the relative difference between the data and the correspond- 
ing theoretical prediction. The hatched band represents the renormalisation- scale 
uncertainty of the QCD calculation. ^ „ 



ZEUS 




Figure 2: Inclusive dijet (a) and trijet (h) cross sections as functions of E^r^^^^y^ 

with the jets ordered in E^^'^^icm- cross sections of the second and third jet were 
scaled for readability. Other details as in the caption to Fig. 1. 
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ZEUS 





Figure 3: The inclusive dijet (a) and trijet (c) cross sections as functions of 

'^LAB '"^^^^ j^^^ ordered in ?7lab-' ''7lab > ''^lab > ''7lab- cross sections of the 
second and third jet were scaled for readability. Figures (b) and (d) show the dijet 

and trijet cross sections as functions of |A?7j^^^j^| between the two jets with highest 

-^THCM- Other details as in the caption in Fig. 1. 
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Figure 4: Diiet cross sections as functions of AEi^^^^j^. 
tions at 0{al) {0{al)) are shown as dashed (solid) lines, 
plots show the relative difference between the data and the 
boundaries for the bins in AE^^^^^'^-^ are given in Table 
caption to Fig. 1. 



The NLOjet calcula- 
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0{a^g) predictions. The 
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Figure 5: Trijet cross sections as functions o/ Ai?^*j^^f^. The measurements 
are compared to NLOjet calculations at 0{al). The boundaries for the bins in 
A_E';|?*j^^j^ are given in TablelSl Other details as in the caption to Fig. 1. 
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Figure 6: Dijet cross sections as functions of . The NLOjet calcu- 

lations at {0{al)) are shown as dashed (solid) lines. The lower parts of the 

plots show the relative difference between the data and the 0{al) predictions. The 
boundaries for the bins in \TipJ.^^(^^\ are given in Table{^ Other details as in the 
caption to Fig. 1. 
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Figure 7: Trijet cross sections as functions of The measurements 

are compared to NLOjet calculations at 0{al). The boundaries for the bins in 
I^PthcmI ^'^^ given in Table\^ Other details as in the caption to Fig. 1. 
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Figure 10: Dijet cross sections as functions o/ | A^'^j^j^^ | . The NLOjet calcu- 
lations at 0{a1) {0{al)) are shown as dashed (solid) lines. The lower parts of the 
plots show the relative difference between the data and the 0{al) predictions. The 
boundaries for the bins in \A(j)^-^(^^\ are given in Tablel^ Other details as in the 
caption to Fig. 1. 
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Figure 11: Trijet cross sections as functions of |A0i^^^^|. The measurements 
are compared to NLOjet calculations at 0{a^). The boundaries for the bins in 
I^'/'hcmI '^'^^ given in TablelSi Other details as in the caption to Fig. 1. 



26 



ZEUS 



lie 

Q. 
X 

CM 



10= 



^ 



10 GeV' < q' < 15 GeV' |Ait>'°^^^^ | ■ 



271/3 




• ZEUS 82 pb'^ 
dijets 

-- NL0iet:0(a2)®c,^_^ 

- NL0iet:0(aJ)®C,3, 

jet energy scalg uncertainty 
□ 1/16<|i2/(Q^+ET)<1 

H 1 1 1 1 1 — I I I 



2x10"' 



10"' 



> 
O 

X 

CM 



10= 



2: 

r 1^ 

^ OE 

3x10"^ 



20 GeV' < q' < 30 GeV' \Aif^^^ I < 27i/3 

• ZEUS 82 pb'^ 
dijets 

-- NLOjet:0(«2)®C,3, 
- NL0jet:0(aJ)®C,3, 




■ jet energy scalg uncertainty 

□ 1/16<h^/(Q^Et)<1 (b) 

— I 1 1 — I — h- 1— I 1 h 



10"' 



2x10"' 



^BJ 



■^rio^ 
> 

O 

Q. 

"xlO^ 

CM 

o 



10 



r 1^ 

^ 



10 GeV' < q' < 15 GeV' |Ait>'°i,'^^ | < 27:/3 




2 _ 



-4 



2x10"' 



10"' 



> 1 0'' F^O GeV < Q' < 30 GeV' lA^i'j;';; I < 271/3- 

O 



x'lO=' 

CM 

o 

■D 

^10^ 



10 

2 

r 1 

^ 

-1 



• ZEUS 82 pb'^ 

trijets 
- NLOjet: 0(aj) ® 




■ jet energy scajg uncertainty 

□ 1/16<h2/(qVEt)<1 (d)- 

— I 1 1 — I — h- 1— I 1 h 



3x10"' 



10"' 



2x10' 



'Bj 



^Bj 



Figure 12: The dijet and trijet cross sections for events with |A0i^^^j^| < 27r/3 
as functions of xbj in two different Q'^-bins. The NLOjet calculations at 0(q!^) 
{0{al)) are shown as dashed (solid) lines. The trijet measurements are compared 
to NLOjet calculations at 0{ag). The lower parts of the plots in (a) and (h) show 
the relative difference between the data and the O(a^) predictions. Other details as 
in the caption to Fig. 1. 
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